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Extracting the Cosmological Constant from the 
Wheeler DeWitt Equation in a Modified Gravity 
Theory 



Remo Garattini 



D | Abstract. We discuss how to extract information about the cosmological constant 

£^ ■ from the Wheeler-DeWitt equation, considered as an eigenvalue of a Sturm-Liouville 

, problem. A generalization to a / (i?) theory is taken under examination. The equation 

is approximated to one loop with the help of a variational approach with Gaussian trial 
wave functionals. We use a zeta function regularization to handle with divergences. 
O , A renormalization procedure is introduced to remove the infinities together with a 

renormalization group equation. 
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1. Introduction 



The Einstein's field equations represent a fundamental set of information regarding the 
laws of space-time. They are represented by 

Rfiv - \g^R + Kg^u = kT^, (i) 

where T^ v is the energy-momentum tensor of some matter fields, k = 8tvG with G 
the Newton's constant and A c is the cosmological constant . The sourceless version of 
Eqs.([H) is simply 

R»v - \g»uR + A c g»v = o. (2) 

It is well known that there exists a huge discrepancy between the observed[lJ and the 
computed value of the cosmological constant. It amounts approximately to a factor 
of 120 orders of magnitude: this is the cosmological constant problem. One possible 
approach to such a problem is given by the Wheeler-DeWitt equation (WDW) [2 J . The 
WDW equation can be extracted from the Einstein's field equations with and without 
matter fields in a very simple way. If we introduce a time-like unit vector u^ 1 such that 
u • u = — 1, then after a little rearrangement, we gel(|| 

M y9 ( 3 d _ 



n = (2k) g^tt v -^rR- 2A C = 0. (3) 



% See Ref.[3] for more details. 
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3 R is the scalar curvature in three dimensions. This is the time-time component of 
Eqs.(T5]). It represents a constraint at the classical level and the invariance under time 
reparametrization. Its quantum counterpart 

H^=0 (4) 

is the WDW equation. If we integrate over the hypersurface E and we define 

ijkl _ V03 ; 



A, 



(2k) G i j k i / !t V/ !T 



2k 



Eq-dlD can be cast into the following form 



V M [gi 



d 3 xA^ 



ri 



(5) 



(6) 



where we have multiplied Eq.fllj) by ty* [gij], we have functionally integrated over the 
three spatial metric and where we have defined the volume of the hypersurface E as 
V = / 2 d 3 Xy/g. Thus one can formally re-write the WDW equation as 

We can gain more information considering a separation of the spatial part of the metric 
into a background term, g^, and a quantum fluctuation, hij, 



lJV[ gij ] tt* [^]/ E d 3 xA s ^[^] 


1 (* 


/ E d 3 xA s 


*> 
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9ij — 9ij 
Thus Eq.(JZJ) becomes 



where A^ represents the i th order of perturbation in hij. By observing that the kinetic 
part of As is quadratic in the momenta, we only need to expand the three-scalar 
curvature J d 3 x^fg 3 R up to the quadratic order and we get 



h 



1 (* 


J E d 3 x 


Ag J + A« + Ag> + . . ." 


*> 


V 









K 



(8) 



(9) 



d 3 X\fg 



-hAh + -h li Ahu - -h ij V l VM+ 
4 4 2 * j 



; + -hRijh ij + jh (rW 



1 

4 J 



(10) 



where h is the trace of and i?*- -* is the three dimensional scalar curvature on- 
shell. Eq.flT]) represents the Sturm-Liouville problem associated with the cosmological 
constant. The related boundary conditions are dictated by the choice of the trial 
wavefunctionals which, in our case are of the Gaussian type. Different types of 
wavefunctionals correspond to different boundary conditions. Extracting the TT tensor 
contribution from Eq.^7} approximated to second order in perturbation of the spatial 
part of the metric into a background term, g^, and a perturbation, h^, we get 



A 



1 



d 3 Xy/g~G 



(2k) K 



-1_L 



, X i X )ijkl + 



(2k) 



;a 2 );^ ( x ,x) 



iakl 



'ID 
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Here G^ kl represents the inverse DeWitt metric and all indices run from one to three. 
The propagator K 1 - (x, x) iakl can be represented as 



x, y 



iakl 



E 



h 



(^) ( y 



2A r 



(12) 



where /i 



(t)± 



x 



are the eigenfunctions of A 2 , whose explicit expression for the massive 



case will be shown in the next section, r denotes a complete set of indices and A (r) are 
a set of variational parameters to be determined by the minimization of Eq. ffTTl) . The 
expectation value of A^ is easily obtained by inserting the form of the propagator into 
Eq. ( TTTT) and minimizing with respect to the variational function A (r) . Thus the total 
one loop energy density for TT tensors becomes 



A 



87TG 



T + 



The above expression makes sense only for uf (r) 
A 2 . 



(13) 

> 0, where Ui are the eigenvalues of 



2. One loop energy Regularization and Renormalization for a / (R) 
theory 



R 



The Spin-two operator for the Schwarzschild metric in the Regge and Wheeler 
representation^, leads to the following system of equations (r = r (x)) 



dx 2 ^ 



. d 2 
dx 2 



+ 



1(1+1) 



1(1+1) 



+ ml (r) h (x) = uhfi (x) 



(14) 



+ m\ (r)\ f 2 (x) = u)l J 2 (x) 

where reduced fields have been used and the proper geodesic distance from the throat 
of the bridge has been considered. Close to the throat, the effective masses are 

ml (r) ~ —ml (M) 

(15) 

ml (r) ~ m\ (M) 

where we have defined a parameter ro > 2MG and tuq (M) = 3MG/tq. The main 
reason for introducing a new parameter resides in the fluctuation of the horizon that 
forbids any kind of approach. It is now possible to explicitly evaluate Eq. (1131) in terms 
of the effective mass. To further proceed we use the W.K.B. method used by 't Hooft in 
the brick wall problem [5] and we count the number of modes with frequency less than 



1, 2. Thus the one loop total energy for TT tensors becomes 



A 



8ttG 



Pi + P2 



167T 1 



E 



mi 



[r)duji 



(16) 



=i J\l m 1( r ) 

where we have included an additional 47r coming from the angular integration. Here, 
we use the zeta function regularization method to compute the energy densities p\ and 
p 2 . Note that this procedure is completely equivalent to the subtraction procedure 
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of the Casimir energy computation where the zero point energy (ZPE) in different 
backgrounds with the same asymptotic properties is involved. To this purpose, we 
introduce the additional mass parameter \x in order to restore the correct dimension for 
the regularized quantities. Such an arbitrary mass scale emerges unavoidably in any 
regularization scheme. One gets 



m\ (r) 



In 



mf (r) 



+ 21n2 



256tt 2 

% = 1,2. The renormalization is performed via the absorption of the divergent part into 
the re-definition of the bare classical constant A 

A ^ A + A div . (18) 

The remaining finite value for the cosmological constant reads 



Ao 
8ttG 



+m 9 




+ 21n2 



( Pl (m) + P2 (fi)) = p T J f (^r).(19) 



The quantity in Eq. ffT9l depends on the arbitrary mass scale /i. It is appropriate to 
use the renormalization group equation to eliminate such a dependence. To this aim, 
we impose that [6] 

1 dA? T (ji) d 



il» 



(20) 



8nG dfx 

Solving it we find that the renormalized constant A should be treated as a running one 
in the sense that it varies provided that the scale is changing 



G 



A (/!, r) = A (/iQ, r) + — (m{ (r) + m\ (r)) 



16vr 

Substituting Eq.(ETJ into Eq.(PSD we find 

Aq(/x ,M) 1_ ' 

8nG 128tt 2 
Eq. fl2^1) has a maximum wherj§| 
1 _ ml(M) 
e Auq 



In 



fM)' 



A (fio 




Gfi 



(21) 



(22) 



(23) 



4/X§ "»vr-»/ 327r 2?re 2- 

The computed cosmological constant appears to depend on the Schwarzschild radius. 
This dependence simply reflects the fact that the chosen background introduce one 
physical scale: the Schwarzschild radius. Nothing prevent us to consider a more general 
situation where the scalar curvature R is replaced by a generic function of R. Therefore, 
we will consider the Sturm-Liouville problem of Eq.(Ej) in the context of a / (R) theor){jj]. 

§ Remark Note that in any case, the maximum of A corresponds to the minimum of the energy 
density. 

|| A recent review on the problem of / (R) theories can be found in Ref.[7]. A more general discussion 
on modified gravities of the type / (R) , / (G) and / (R, G) where G is the Gauss-Bonnet invariant, can 
be found in Rcf.|8|. 
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3. One loop energy Regularization and Renormalization for a generic / (R) 
theory in a Hamiltonian formulation 

In this section, we report the main steps discussed in Ref . [9] for a / (R) theory in 
connection with the Sturm-Liouville problem of Eq.([7j). Although a / (R) theory does 
not need a cosmological constant, rather it should explain it, we shall consider the 
following Lagrangian density describing a generic f(R) theory of gravity 

C = y/=jj(f(R)-2A), withf^O, (24) 

where / (R) is an arbitrary smooth function of the scalar curvature and primes denote 
differentiation with respect to the scalar curvature. A cosmological term is added also 
in this case for the sake of generality, because in any case, Eq. (|24l) represents the most 
general lagrangian to examine. Obviously /" = corresponds to GR. The generalized 
Hamiltonian density for the / (R) theory assumes the fornjjjl 

1 r -y/gf (R) (®R - 2A C - SKijlCi + K 2 



n = — 

2k 



where 



and 



V 



+V(V)+2g^(^f(R)) l ..-2^K ij 
- QVgf (R) 



(25) 
(26) 



V(P) = y/g[Rf{R)-f{R)\. (27) 

Henceforth, the superscript 3 indicating the spatial part of the metric will be omitted 
on the metric itself. When / (R) = R, V(V) = as it should be. Eq.( l25l) becomes 



n = f (R) 
1 

+ 2k 



(2k) G ijki TX %3 TX 

yfgfiR) [2K l3 K 



V(V)+2gV( x ^f'(R)) l ..-2pVK l 



Since 



(28) 
(29) 



then we obtain 

n = f(R) 
i 



(2k) G ijki Tx l H 



2k V 



2A, 



+ 2k 



2^K VJ K^ (f (R) - 1) + V(V) + 2g« (rff (R)) {ij 
and transforming into canonical momenta, one gets 

ij^ki Vd 



(30) 



n = f(R) 



(2k) G ijkl rfH k 



y* m R 
2k V 



2A, 



% See Ref.fTu] for technical details. 
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(f (R) - l)+i- [V(V) + 2gv (y/gf .(31) 



By imposing the Hamiltonian constraint, we obtain 

f'(R) 
l 



(2k) G ijkl ^-K kl - 



+ 



2k 



V(V) + 2g»(y/gf (R) ) 



■2(2k 

-f (R) Vg 



ii 7T 

Gnki^n + — 



A, 



(/'(^)-l) 

(32) 



If we assume that /' (R) ^ the previous expression becomes 



2k 



+ (2k) 



7T 

GijklTT 13 ^ + — 



+ 



2Kf'(R) 



-A r 



Hf'(R)-i) 



f'(R) 



Now, we integrate over the hypersurface E to obtain 
d'xl (2K)G, ilH ^7T kl -^R 



+ / d 3 x 



The term 



s 2«/'(i2) 



1 



+ (2k) 
V(V)+2gV(^gf(R)) l . j 



7T 

GijklTT lj TT + — 



(33) 
2(f (fl)-l) 



A, 



- ' d3x pm9 lj (Vgf (R)) {ij 



d 3 x^g-. (35) 



(36) 



appears to be a three-divergence and therefore will not contribute to the computation 
The remaining equation simplifies into 



drx 



(2k) G^tt" - ^R 



V(V) 



K 

K 



2k 

J d 3 x^g. 



+ (2k) 



GijklTC tj TC H + 



7T 



Hf'(R)-i) 



f'(R) 

2k f (R) J k is^ v ^ ^ 

By a canonical procedure of quantization, we want to obtain the vacuum state of 
a generic f(R) theory. By repeating the same procedure for the generalized WDW 
equation Eq. (j3~T|) . we obtain 



1 (* 


kd 3 x 




*> 


V 







*) 2k 2 (/' (R) - 1) (* Is rf3;E G ijkl ^TT kl + 7T 2 /4 



+ 



1/ f(i?) 



1 (fr|/ s cfe^)/(2^f QB))|fr) _ Ac 

V (#|#) ~ K 



(*|*> 



(3f 



From Eq. (138jl . we can define a "modified" operator which includes /' (i?). Thus, we 
obtain 



kd 3 x 



A (2) 



«7 , «(/' (#)-l)(*|/ s d 3 o;[7r 2 ]i*> 
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Ag ) /(fl) = (2K)/i(i2)G yfc ,7r^7r H 



Ac 



2k 



ft, = 1 + 



2 [f (g) ~ 1] 

f'(R) 
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(39) 

(40) 
(41) 



and where 3 R hn is the linearized scalar curvature whose expression is shown in square 
brackets of Eq. flTO]) . Note that when f (R) = R, consistently it is h(R) = 1. From 
Eg. (1591 . we redefine A c 



a: = a c + 



1 f j3 -Rf'(R)-f(R) ,. n . 
2V {m : X -2V.L (l1 ^ Ar) • (42) 

where we have explicitly used the definition of V(V). In the same spirit of the previous 
section, we restric t th e analysis to the contribution of physical degrees of freedom, 
namely TT tensors . Thus Eq. ( 122|) becomes 



A(,(/x ,r-) m 4 (M) 



8ttG 



128tt 2 



In 



'ml (MY 



+ 



(43) 



Now, we compute the maximum of A , by setting x = m\ (M) /4/i 2 ,. Thus A becomes 



Aq (a*o, x) 



fln(x) + i 



TT 



As a function of a;, A (fio, x) vanishes for x 
x G 



and x 



exp 



(44) 



-|) and when 



0, exp |) , Aq(^O)^) > 0. It has a maximum for x — 1/e equivalent to 

(45) 



wIq (M) = 4/ig/e and its value is 
A (A*o,a? " 



2vre 2 



or 



1 /" J3 r- R f ( R ) ~ f ( R ) 



'h (R) 

Isolating A (fMo, x), we get 

A (/i , x) = yjh {R)^ ~ ^7 / s 
Note that Ao (no, x) can be set to zero when 

h[K) 2^~ 2Vh d XV9 f'(R) 



Gl4 
2vre 2 ' 



Rf (g) - / (R) 



(46) 



(47) 



(4c 



+ For a complete derivation of the effective action for a / (R) theory, see Ref. [TT] . 
* By a canonical decomposition of the gauge part £ Q into a transverse part £j with V a £j = and a 
longitudinal part £a with ^Jj = V a i/>, it is possible to show that most of the contribution comes from 
the longitudinal part (scalar). Evidence against scalar perturbation contribution in a Schwarzschild 
background has been discussed in Ref.[12j. 
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Let us see what happens when / (R) = exp (—aR). This choice is simply suggested by 
the regularity of the function at every scale. In this case, Eq. (1481) becomes 



3aexp(-afi) + 2 g^ = 1 , 
\ aexp(— an) ire^ aV Jt, 

For Schwarzschild, it is R — 0, and by setting a = G, we have the relation 

^ = 7T\/(3G + 2)G- (50) 

Like the case of / (R) = R, the purpose of this calculation is related to Eq.©, which 
can be applied to different backgrounds, i.e. de Sitter, Schwarzschild-de Sitter, etc. The 
result of this process will be a spectrum built in terms of vacuum fluctuations. Once 
the " Ground State" of this spectrum will be identified, then there will be a chance to 
have an approach for explaining the so-called Dark Energy. Note that in this approach, 
there is no evolution in time, so some pathologies of the type "Big Rip" do not come 
into play at this stage. 
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